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We consider multidimensional cosmologies in even-dimensional space-times [D — 2n) containg perfect fluid and a 
multidimensional generalization of the Maxwell field Fai - a„ preserving its conformal invariance. Among models 
with an isotropic physical 3-space some integrable cases are found: vacuum models (which are integrable in the 
general case) and some perfect fluid models with barotropic equations of state. All of them contain a component 
of the F field appearing as an additional scalar in 4 dimensions. A two-parameter family of spatially fiat models 
and four one-parameter families, including non-spatially flat models, have been obtained (where the parameters are 
constants from the fluid equation of state) . All these integrable models admit the inclusion of a massless scalar field 
or an additional fiuid with the maximally stiff equation of state. Basic properties of vacuum models in the physical 
conformal frame are outlined. 



1. Introduction 

This paper continues a study of exact solutions to the 
gravitational field equations including a generalization 
of the Maxwell field to 2n -dimensional space-times 
(n = 2,3,...), began in Ref. 11. UnHke the conven- 
tional straightforward use of the Maxwell field in D di- 
mensions, this generalization preserves the invariance 
of the field under multidimensional conformal transfor- 
mations. 

This generalization (Generalized Maxwell Field, 
GMF) is represented by a (_D/2)-form, where D = 2n 
is the space-time dimension: 



F = dU, or Fa,. 



nAilJ A2...A„] (1) 

where J7 is a potential (n— l)-form and square brack- 
ets denote alternation. The field F is invariant with 
respect to the gauge transformation 



U 



U + dZ 



(2) 



where Z is an arbitrary {n — 2) -form; in 4 dimensions 
this is the conventional gradient transformation of the 
Maxwell field. 

The GMF was introduced in Refs. ||, |], where 
some particular cosmological models with the GMF 
were studied; in ^ we discussed in some details sev- 
eral families of spherically symmetric solutions, with 
combinations of electric, magnetic and quasiscalar type 
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components of the GMF, with possible inclusion of a 
massless, minimally coupled scalar field (p = . 
Some cases of interest were found, namely, black holes 
with a somewhat unusual mass and charge dependence 
of the Hawking temperature. 

This paper presents some classes of exact cosmo- 
logical solutions, including those of ^, ^] as special 
cases. Just as in those papers, only Friedmann-like 
models, isotropic in the three physical dimensions, are 
considered, therefore the GMF can have only "quasis- 
calar" components, behaving as scalars in the physical 
space-time. The solutions found here contain only one 
such component. Further generalizations with a num- 
ber of such components and even a number of p-forms, 
as discussed in field models like M-theory (see e.g. jH) 
and references therein), are possible. 

The aim of this paper is basically to indicate the 
integrable cases. A more detailed description of the 
solutions will be given elsewhere. 

Only the behaviour of vacuum solutions is briefly 
outlined, some of them exhibiting attractive features. 
In particular, solutions with a hyperbolic 3-space, have 
a linearly expanding physical space at their late-time 
asymptotic, with internal-space scale factors tending to 
finite constant values which (bearing in mind the usual 
assumption that the internal spaces are compact) may 
be arbitrarily small. 

Throughout the paper capital Latin indices range 
over all D space-time dimensions, Greek ones take the 
values 0,1,2,3, indices like ai, . . . ,an refer to internal 
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spaces and i,j enumerate the factor spaces. 

2. Generalized Maxwell field and 
3-isotropic cosmology. 
Field equations 

Consider general relativity in a Riemannian space-time 
(D = 2n, n > 2), in the presence of the GMF and 
perfect fluid (matter), with the total action 

F^ = F^'-^-Fa,...a^. (3) 

where R — ^ R is the scalar curvature corresponding 
to the £) -metric gAB , 9 = Idetgy^sl, k is the D- 
dimensional gravitational constant, F is the GMF and 
Lm is the matter Lagrangian. 

The -sector of the action is invariant under con- 
formal mappings of {qab ^ f{x^)gAB) provided 
the components i^Ai...yi„ are unchanged in such trans- 
formations (the potential U from (^ may experience a 
gauge transformation (^) where Z is some (n— 2) -form 
depending on the conformal factor / ) . 

The field equations are 



Ga = Ra ^ \^aR — ^'^T^A 1 

y^pAA,...A„ = 0, 



(4) 
(5) 



where the energy-momentum tensor (EMT) is a 
sum of contributions from F and matter. The EMT 
of the F field is 

Tj = (-l)"-^(nF^c....c„^^^^-^" - \5lF^). (6) 



An alternative form of (0) is 



D-2 



B rjiC 



(7) 



where is the EMT of aU matter except the F field. 

m 

In what follows we will assume that the physical 
3-space is isotropic and the GMF has only one (up to 
index permutations) nontrivial component compatible 
with such isotropy, namely. 



F 



(8) 



where the index corresponds to time and all the oth- 
ers to extra dimensions; moreover, we restrict ourselves 
to the simplest space-time structure compatible with 
such an F field. Evidently, among the 2n — 4 internal 
dimensions Eq. (^ distinguishes two subsets, with and 
without the F field, or, in other words, n — 1 coor- 
dinates from the remaining n — 3. Thus for n > 3 
there must be at least two internal factor spaces with 
different dynamics. 



Accordingly, we consider V with the structure 



V^^ = IR X Ml X M2 X M3, 



(9) 



where IR corresponds to time, Mi is the external 
isotropic 3-space (of curvature K — 0,±1), M2 and 
Afa are Ricci-flat internal spaces, dimM2 — n — 1, 
dim Afa = n — 3 , so that M2 is parametrized by the 
coordinates whose numbers belong to in (H) and 
Afa to the remaining ones. 

The metric is assumed in the form 



s'f'^dsl 



s'f^dsl 



(10) 



where 7, Pi are functions of the time variable r and 
dsi are the r -independent metrics in Mi. 

For the perfect fluid we assume a pressure isotropic 
in each factor space, so that its EMT is 



Ta = diag(/9, [-pi]3, [~P2]n-l [~P3]n- 



(11) 



where p is the density, pi are the pressures in the cor- 
responding directions and the symbol [x]k means that 
a quantity x occupies k positions along the diagonal. 
We admit that p and pi are connected by the equa- 
tions of state 



Pi = [a-i - l)p, 



Qi = const. 



(12) 



so that a physical range of is between (vacuum- 
like state) and 2 (stiff matter) . 

It should be noted that the 4-dimensional part of 
the _D -metric ( p^ ) does not coincide with the metric to 
be used to interpret actual cosmological observations. 
Indeed, real space-time measurements (such as Solar- 
system experiments or cosmological redshift measure- 
ments) rest on the constancy of atomic quantities (the 
atomic system of measurements). Thus, the modern 
definition of reference length is connected with a cer- 
tain spectral line, determined essentially by the Ryd- 
berg constant and ultimately by the electron and nu- 
cleon masses. Therefore it is reasonable to apply as a 
physical one such a metric (to be denoted .g*^, ) con- 
formally related to 5^,^ with which masses of bodies of 
nongravitational matter, like atomic particles, do not 
change from point to point. This happens if the mat- 
ter Lagrangian appears in the total Lagrangian in 4 
dimensions without any x'' -dependent factor. 

The choice of 5*^ depends on how Lm appears in 
the original action. In our case (||), just as in Ref. it 
is easy to check that the metric in the "atomic gauge" 
should have the form 

= e-/'ffM- (13) 

where e'^ is the volume factor of extra dimensions, in 
our case a = {n—l)f32 + {n — 3)P3. Thus, in particular, 
the physical 3-dimensional scale factor a(r) is 

air) = exp [/3i + l{n - 1)^2 + lin - S)^^]. (14) 
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We are now ready to write down the field equations 
explicitly. Let us use the harmonic time coordinate t , 
so that 



7 



(r) =^^7V,/?.(r), 



= dim Mi. 



(15) 



The Maxwell-like equation (||) for the only nonzero 
GMF component -Fb4...n+2 is easily integrated giving 

- = q/V9^ 1 = const (16) 

where due to (|l|) ^5 = e^'' ■ So the EMT of the F 
field may be written as 



F 



xdiag(l,[-l]3,[l]„_i,[-l]„_ 



(17) 



For perfect fluid the conservation law "S/ bT^ = 



after integration leads to 

p = Pa cxp f - ^ Niai|3^ j . 

^ i=i ' 



po = const 



(18) 



and its EMT is presented in the form 

111 

X diag(l, [1 - ai]3, [1 - a2]„-i, [1 - az\n-z) 

(19) 

where we have introduced the function 

3 

2y - ^(2 - a,)NA. (20) 
1=1 

With ( p^ ) and (p^), the nonzero components of the 
Ricci tensor are 



J- . ^ITl ■ 



e2'5'X,(A^,-l) 



(21) 



(22) 



where if^ — 0,±1 are the factor space curvatures; we 
will assume K2 = = and only for the physi- 
cal space treat all the three variants, corresponding to 
different Friedmann-Robertson- Walker models. 

The spatial components of the Einstein equations 
now read: 



2X1 e 



P2 = -g'c2-+A2Kpoe'^ 
Ps = 



(23) 
(24) 
(25) 



where 

X = {n- l)/32, 



/ii — a,; 



A^2 
2n- 2' 



Q2 = ^nlnq^. 



(26) 



Eqs. (P^)-(p5|) have a first integral (the energy in- 
tegral), coinciding with the (^) component of the Ein- 
stein equations: 



kpq e 



2y 



(27) 



In what follows we try to solve Eqs. (|2|)-(|5|); ^) 
is used to obtain an additional relation between the 
integration constants and simultaneously to verify the 
correctness of the solutions. 

3. Vacuum solutions with the F field 

p = Pi = one easily solves 



In the case po = 

Eq. my. 



'("-1)^^ = ^^/{n-l)Q^ cosh[/i(7 



T2 



)] (28) 



with the integration constants h > and T2 . As now 
P2 + /^a = , (33 is also easily found: 



P3 



-/32 +6T + 63, 



6, b. 



const. 



(29) 



Lastly, (23) can be combined with ( p4| ) and ( |2^ ) to give 
(7-/3i)---)-4ifie2^-2'3i ^0, (30) 
whence 



p/3i-7 



2s(fc,r), 



/^i -1 (fc e R); 
==0, (fee R); 



2fc~icoshfcT, = 1, (fc>0), 

(31) 



where 



s{k, t) 



k ^ sinh fcr, fc > 
T, fc = 

fc~^ sin fcr, fc < 



(32) 



fc — const and one more constant is dominated by 
shifting the origin of r . The integration is completed. 
A substitution to the energy integral (27) gives: 

-h^ + (n- l)(n- 3)6^ 



3/1: = 



K 



dcf 



n-1 
2 



fc signfc, K = -1; 
fc2, i4: = 0,+l. 



(33) 

1 we have fc > and 



so that, in particular, for Ki 
e"^ = (2/fc)sinhfcr. 

The physical components of the metric are (see 
®, ©)• 

e^^ = exp [3z — P2 — \ {n — 3)6r]; 



5oo 

2(r) = exp[z-/32-i(n-3)6T]. 



(34) 



The model properties strongly depend on the di- 
mension 2n and specific values of the integration con- 
stants. Thus, the cosmic proper time t — J exp(7*)dT 
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is finite or infinite depending on the convergence or 
divergence of the integral, etc. One can conclude the 
following. 

1. The harmonic time t e IR for Ki = 0, +1 and 
(with no generality loss) r e (0,oo) for Ki — —1. 

2. The internal-space scale factor e^^ tends to zero as 
T —^ ±oo , i.e. at both ends of the evolution process 
for Ki = 0, +1 and at one end for Ki = — 1. In the 
latter case, P2 tends to a finite limit as r ^ . 

3. Another scale factor e/^a (existing only for n > 3) 
tends to infinity either as t —^ +00 , or as t ^ —00 , 
or in both limits (if |6| < h/2). Again, it is finite as 
T^Q {Ki = -1). 

4. For n = 3 {D = 6), Eq. (||) gives just 3P = 
2h^ , therefore the solution behaviour is uniquely char- 
acterized in terms of the proper time t. Thus, for 
Ki — (spatially flat world) the expanding model has 
the aysmptotics (where arbitrary constants are chosen 
with no loss of generality) 



„(-V3/2+l/4)/i|r| 



t 



a{t) ^ e(-Vi/6+i/4;^|^| 



initial singularity, and 



0, 



i0.162 



(35) 



a{t) ^ e(Vi/6+i/4)/.|r| 



.0.447 



1/2 



,-0.34 



(36) 



— expansion in physical space and contraction in extra 
dimensions. 

For Ki — +1 (spherical world), the time symmetry 
of the evolution is broken only by the shift T2 in e^ 
which does not affect the asymptotics. The evolution 
proceeds between two singularities like (35). 

For Ki = — 1 (hyperbolic world) the model has a 
singularity like (Ba) at r = cxd ; on the other hand. 



t 



ait) 



00; 

t as 



0, 



(37) 



i.e. there is a linear expansion (contraction) asymp- 
totic with the internal scale factor e''^ tending to a 
finite limit. 

5. For n > 3, the advent of one more constant b makes 
the evolution of a(t) more diverse. One can observe, 
in particular, that, for any set of constants, there is 
always an asymptotic like (|37| ) for Ki — —1 and at 
least one singularity like (|35D (with, in general, other 
numerical characteristics) for K — +1. 

It should be stressed that the behaviour of a{t) 
and even the finitcness or infinitcncss of time intervals 



crucially depend on the choice of the conformal frame. 
Thus, the present vacuum solutions for n — 3 coincide 
with those of Ref. [Q , but the conclusions are different 
since there the authors considered the metric (7^,^ as 
the physical one. 

In what follows we only indicate the integrable 
cases, postponing their detailed analysis to further 
publications. 

4. Solutions with matter: Ki = 

The unknown function y is excluded from Eq. (|2j) in 
two cases: 

1. if y is proportional to x (i.e., to /32), 

2. if yl2 0. 

1. This possibility is realized if ai = 03 = 2 (stiff 
matter equation of state in Mi and M3 . Then we 
have 



/3i = /33 = -P2 



(38) 



so that it is sufficient to integrate (M for (32 , then /3i 
and will be easily found. Eq. ( p4[) reduces to the 
form 



/32 + g'e2("-i)'3^ + (l-iai) 



Thus it belongs to the type 



-,2a; 



B2 e^""" = 0, 



Bi, B2, c — const. 
Its first integral is 

i;2 + Bi e^^ + {B2/C) e^"'' = const 



(39) 



(40) 



(41) 



which is evidently integrable by quadratures. The so- 
lution is expressed in elementary functions in the cases 
c = 1/2, 1, 2. By (39) this corresponds to 



0,2 



1, 0, -2, 



(42) 



the first two values belonging to the physical range. 

So we have obtained a class of solutions which can 
be labelled 



(2n,0,F| 2,a2,2), 



(43) 



where the first three symbols designate the dimension, 
spatial curvature (in Mi ) and a nonzero F field, while 
the remaining three indicate the coefficients . 

2. The relation ^2 = is explicitly written as 

2 + (n- 1)02 = 3ai + (n- 3)03. (44) 



So the equation for /32 is integrated in the form (28). 
Then, if we combine Eqs. (|2^)-(^) to form y in the 
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left-hand side, it turns out that, in the right-hand side, 
the coefficient before e^^ vanishes and we obtain 

2y = [3(2-ai)(ai-a2) + (n-3)(2-a3)(a3-a2)] e^^ 

(45) 



which is easily integrated like (30). So we know P2 
and y . Lastly, there is a linear combination of (23)- 
( p5[ ) whose right-hand side is zero, namely, 

(03 - a2)/3i + (aa - ai)/32 + (02 - ai)/33 0. (46) 

In the general case Eqs. ( |45|) a nd ( ^6|) are independent, 
so after trivially solving ([46|) the integration is com- 
pleted. In the exceptional case when they are not, one 
can find /3i from ( p3| ) , since both x and y are already 
known. 

Thus for Ki = we have obtained two families of 
solutions: (|4^ ) (one-parameter) and 



(2n, 0,F I ai,a2,a3), 



(47) 



constrained by (^4|) (two-parameter). Their only com- 
mon element is the trivial case of all Ui — 2, when all 
Ai = and the solution is reduced to vacuum up to 
the relation between the integration constants that fol- 
lows from the energy integral. This case of stiff matter 
in the whole space coincides with that of a minimally 
coupled scalar field (^(r) , to be discussed later. 

5. Solutions with matter: Ki = ±1 



The same combination of (|2^)-(^) that has led to 
Eq. (pfl) in the vacuum case, for po 7^ reads: 



2y 



(48) 



One can again select two cases of its integrability: 

1. ai = 2 [then (H) coincides with (js^)] and 

2. y is proportional to (7 — /3i) and ( ^ reduces to 

1. ai = 2, stiff matter in physical space. We are 
left with Eqs. ( p^ and (^ ) with two unknowns x and 
y. Eq. ( p^ is easily integrated (reduces to (28)) when 
A2 — 0, that is 



{n — l)a2 — (n — 3)03 = 4. 



(49) 



This condition is of interest only for n > 3, since for 
n = 3 it reduces to a2 = 2, i.e., the trivial case ai = 
a2 = 2. 

For n > 3 , the function y can be found from an 
equation similar to (|4^) but with another constant be- 
fore . This completes the integration. 

We thus obtain the 1-parameter family of models 



(2n,±l,F I 2,02,03), 71 > 3 



(50) 



constrained by (^491). 

For oi = 2, n — i, there is no 13^, and Eq. ([2^ ) 
takes the form 



2/32 



-2g2e4^^+Ajpo(a2-2)e2^'^(2-a,) 



(51) 



belonging to the type (^). It is integrable in elemen- 
tary functions for 02 = 1 (no pressure) and 02 = 
(vacuum- like state). Thus we obtain the models 

(6,±l,F|2,a2). (52) 

Another case with an equation like (|ll) exists for 
n > 3: if oi = 03 = 2, then 

2y= (n-l)(2-a2)/32. (53) 

Thus Eq. (|^ is solved; then, knowing the combination 
7 — /3i and (32 (hence also y ) , it is straightforward to 
find /33 from (^5|). Thus we obtain a family of models 
similar to (E3|) 



(2n,±l,F I 2,02,2) 



n > 3. 



(54) 



2. oi 7^ 2; to solve (p^), let us require that the two 
3- vectors of coefficients by j3i in the expressions for y 
and 7 — /3i be parallel: 

3(2 -oi) \ 
(n-l)(2-a2) ^ I n-1 I , (55) 
,(n-3)(2-a3)/ 
whence 



02 = 03 = |oi — 1. 



(56) 



In this case also automatically A2 — — Q, so that 
(32 and fi^ are found quite easily, thus completing the 
integration. We get the models 



(2n,±l,i^ I 01, foi - 1, §01 1). 



(57) 



Thus we have obtained 4 one-parameter families of 
models with nonzero spatial curvature. In each case, 
to obtain a final form of the solution, it is necessary 
to substitute the results of integration to the first in- 
tegral ( p7| ) to find a relation between the integration 
constants. 

6. Inclusion of a massless scalar field 

So far we have been working with the action (^ . How- 
ever, all the above solutions can easily include a mass- 
less, minimally coupled scalar field tp = (^™'" with the 
action 



(58) 



Then, in the space-time considered, the scalar field 
equation for (p = (^(r) is immediately integrated to 
give 



ip = ipo + (Pit; 



PQtPi — const. 



(59) 
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The EMT of 0""" coincides up to a constant factor 
with that of stiff matter (when aU = 2). 

It is easily verified that 0™'" does not contribute 
to the right-hand side of Eqs. and therefore in no 
way affects the solution process described above. It 
only affects the system through its contribution to the 
first integral where one should add the term K(f1 
to the right-hand side. The same applies if one adds 
one more component of matter in the form of a perfect 
fluid with stiff equation of state (ai = 2): then the 
addition is just Kpo2 = const . This addition certainly 
changes some properties of the solutions. 

Due to the conformal invariance of the F field, D- 
dimensional conformal transformations ||^, |j lead to 
models of generalized scalar-tensor theories of gravity. 
Vacuum solutions are then generalized in a straightfor- 
ward way, whereas perfect fluids will in general acquire 
a nontrivial interaction with scalar fields. 
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